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Abstract 

The vector of weights of an interpolatory quadrature rule with n preas- 
signed nodes is shown to be the least-squares solution uj of an overde- 
termined linear system here called the fundamental system of the rule. 
It is established the relation between w and the minimax solution z of 
the fundamental system, and shown the constancy of the oo-norms of 
the respective residual vectors which are equal to the principal moment 
of the rule. Associated to to and z we define several parameters, such 
as the angle of a rule, in order to assess the main properties of a rule 
or to compare distinct rules. These parameters are tested for some 
Newton-Cotes, Fejer, Clenshaw-Curtis and Gauss-Legendre rules. 

Key-words: Quadrature rule, least-squares solution, minimax solution, prin- 
cipal moment. 
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1 Introduction 

We establish several connections between the computation of the weights u) 
of an interpolatory quadrature rule with n given nodes and the least-squares 
and minimax solutions of a certain inconsistent linear system of (n + 1) 
equations in n unknowns. This system, hereafter called the fundamental 
system of the rule, is obtained by applying the undetermined coefficients 
method to a convenient basis B of polynomials, which depends on the given 
abscissas. We prove that the vector of weights uj has the property of having a 
constant residual norm, which is independent of the norm used, equal to the 
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so called principal moment of the rule. General references for interpolatory 
quadrature rules and its applications are, for instance, [8], [T7], [U], [9], [18]. 

In order to approximate the integral /(/) = f(x) w(x)dx (with w(x) a 
nonnegative weight function, continuous and integrable on the open interval 
(a, b)) we assume that an ordered set M of nodes is given, say t\ < £2 < 
. . . < t n (with ti belonging or not to the interval of integration). We consider 
the quadrature rule 

n 

Q„(/) = (1) 

i=l 

The entries of the vector uj = [ui, . . . ,oj n ] T are the unknown weights of the 
rule. The rule (pQ) is assumed to be interpolatory, meaning that it is exact 
at least for all polynomials p of degree less or equal to n — 1: 

Qn{p) = I{p), Vp G P n -l, 

where Vk denotes the vector space of polynomials of degree less or equal to 
k. The rule Q n (f) is said to have degree of exactness d (or simply degree d) 
if it is exact for all polynomials of degree < d (i.e., for all p 6 Vd) but it is 
not exact for some polynomial in Vd+i- 

It is well known that any interpolatory rule Q n (f) has at least degree n — 1 
and at most 2n — 1 (see for instance [IT], Ch. 7). Consequently, the respec- 
tive weights lo can be obtained using the so called undetermined coefficients 
method applied to a basis of V2 n-i- Traditionally one considers the standard 
basis 4>i(t) = t l , with i = 0, 1, . . . , (2n — 1), and the undetermined coeffi- 
cients method applied to this basis leads to the following overdetermined 
linear system of 2 n equations in n unknowns 

Qn{(t>o) = Mo 
Qn{4>l) = Ml 

■ • (2) 

Qn{4>2n-l) = M2n-1, 

where the moments \X{ are assumed to be finite, and equal to 

w = J(0 i (t))= [ faifyw^dt, i = 0,l,...,(2n-l). 

J a 

Since the degree of the rule is generally unknown, one cannot decide a priori 
what is the relevant set of n + 1 equations in the system ([2]) (see [6], Ch. 2, 
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for the discussion on the choice of relevant equations from an overdetermined 
system). Furthermore, the matrix of the system is generally dense and so 
is natural to ask whether there exists another basis in which the system (|2|) 
has a simpler form. In our work [15] we show that choosing a basis B as in 
Definition 12. 1\ we easily replace ^) by a relevant linear system with n + 1 
equations in n unknowns. Such system will be called here the fundamental 
system of the rule, since it is univocally determined by the set of nodes j\f 
which in turn determines the basis B. Thus, the choice of this particular 
basis B may be seen as the right way to extract the relevant set of equations 
from the system ([2]). Otherwise (taking the standard basis of polynomials) 
the computation of a "solution" of ([2]) relies on heavy tools such as Polya's 
algorithm [19], detailed in [6] (see also [1]). 

In Section[2]we study the main properties of the fundamental system F uj = c 
associated to the rule Q n (f)- This system has the remarkable property of 
containing a triangular subsystem Aui = c, from which the weights w can be 
computed by backward substitution. The last entry of the vector c contains 
a non null quantity /j,q which we call the principal moment of the rule. The 
computation of the parameter \iq gives automatically the degree of exactness 
of the rule (see the work of the author |15j). Therefore, the principal moment 
can be seen as a kind of rule's signature. 

The main results are given in Propositions 12. ll and 12.21 The first proposition 
shows that the least-squares solution of the fundamental system coincides 
with the vector of weights uj, and in Proposition 12.21 (a) we prove that the 
minimax solution z of the fundamental system differs from the least-squares 
solution a; by a correction vector r. This correction r (and so the minimax 
solution) is recursively computed by solving a certain triangular system. 

Several parameters related to the least-squares ui and minimax z solutions 
are defined in order to preview the asymptotic behavior of a rule or to 

compare distinct rules. For a fixed number n of nodes, the first two main 

ii ii*M 
quantities to be considered are respectively N u = \\uj\\i and N* = || z \\i. 

Other interesting parameter is the angle between ui and z, which we call the 
angle of the rule. 

Let us now explain the heuristic behind the choice of these parameters. The 
most remarkable property is that uj and z have residual vectors r(uj) and 
r(z) of constant norms (see Proposition 12.21 (b)) in the following sense 

IkMllp = ll r (^)lloo = \fiQ\, for anyp = oo,l,2,3,... . 
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This property allows to preview that a convergent rule should have their 
parameters N u and N* aproaching zero as n goes to oo. Moreover, the 
angle of a rule should also be close to zero for n sufficiently large. As 
observed in Section [3l the behavior of the norm parameters N u and ./V* for 
the Newton-Cotes rules do not decrease as n increases (see Figure [Q, and 
so one concludes that the Newton-Cotes rules cannot be convergent. 

In Section [3] we test and compare the behavior of several parameters 
for the Fejer (F), Clenshaw-Curtis (CC) and Gauss-Legendre rules (GL), 
for n varying from 2 to n = 17. We give numerical and graphical evidence 
that the Fejer and Clenshaw-Curtis rules have an angle close to the angle 
of Gauss-Legendre rule. This is a plausible explanation for a phenomenon 
observed by several authors, in particular by Lloyd N. Trefethen in [20]: in 
spite of its low degree of exactness, for large n, the CC rule is almost as 
accurate as GL. 

The least-squares/minimax approach to interpolatory quadrature rules gives 
a great deal of information to assess the main properties of a rule, with no 
explicit need to appeal to interpolation theory. Indeed one can deduce an 
error expression, E n (f) (see Proposition 12.11 (c)). in the case the integrand 
function is assumed to be sufficiently smooth. Moreover we do not need to 
make any assumptions on the smoothness of the integrand function / (apart 
from the existence of 1(f)) to compute the principal moment of a rule as 
well as other relevant parameters such as the residual norms or the angle of 
the rule. 

In conclusion, the minimax/least-squares approach to interpolatory quadra- 
ture here discussed gives rise to interesting analytic and geometric aspects 
opening new research directions. In particular, we think it is worth to further 
explore the role of the referred parameters on the fundamental questions of 
the convergence and accuracy of an interpolatory rule. Another direction of 
research is to relate the least-squares/minimax approach with the available 
methods for simultaneous computation of nodes and weights using Jacobi 
matrices (see (T3] and [14], Ch. 10). 

2 The fundamental system of a rule 

Given an ordered set M of nodes we construct a polynomial basis B = 
Bj\f, defined in ([3]) and (j4]), dependent on the nodes. The undetermined 
coefficients method is then applied to B leading to an overdetermined system 
T = J-b,Ni with n + 1 equations in n unknowns. 
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Definition 2.1. (Canonical basis) 

Given a set N of ordered nodes t\ < £2 < • • • < t n , the polynomials 

MX) = 1 ,gN 

4>j = (fij-i(x) x (x — fj), 1 < j < n — 1 

(where the polynomial <f>j(x) has degree j) is a basis of V n —\.This basis is 
here called the canonical basis associated to N '. 

To form a basis of V n +k (where k will be specified in accordance with the 
particular rule under consideration), we complete the canonical basis with 
the polynomials q n (x), q n+ \(x), . . ., q n +k(%) defined by 

q n (x) = <t> n -l(x) (x - i n ) / 4 x 
qj(x) = q,j-i{x) (x — t r ), j > n + 1, with r = j (mod n) 

Remark 2.1. Note that for any j > 1, the zeros of the polynomial 4>j(x) are 
exactly t\ , t% , . . . , tj . Also the nodes t\ , t% , . . . , t n are zeros of each polynomial 
qi{x) given in Q, for any i > n. 

The canonical basis and the extended basis have been used in a previous work 
for the simultaneous computation of the degree of a rule and its weights via 
the undetermined coefficients method. We recall now the main result in |15| . 

Theorem 2.1. fThJ 

Let M be the set of given nodes, t\ < t-i < . . . < t n , and consider the basis 
Bo = {<j)o(t),. . ., <p n -i(t)} ofV n -i, and the basis B = B_\f = B Q U {q n (t) , 
q n +k(f)} ofV n+ k, where (pi and qi are defined in ([3]) and @. Denote 
by Hi the moments fa = <fii(t) w(t) dt, i = 0, . . . ,n + k. 

(a) The undetermined coefficient method applied to the basis Bo deter- 
mines (uniquely) the weights uii of the rule ([T]) for approximating the 
integral 1(f) = J f(x)w(x)dx. These weights are 

Wn = - 



4>n-l{tn) 

ft-i-EL+i^-ifek . 101 (,) 

u)i = — — , i = n-l,n-2,...,l. 

<Pi-i{ti) 

(b) The undetermined coefficient method applied to the basis B, deter- 
mines the degree of exactness d = deg(Q n (/)) as being the integer 
d > n — 1 for which 

Hj = I(qj) = 0, for all n< j <d and n d+l = I(q<i+i) / 0. 



5 



(c) If f is of class C d+1 ([a, b]), the error expression of the rule is 
E Q M) = c n f( d+1 XO, with c- W+1 

for some £ G (a, b). 



{d+l)V 



(6) 



For any polynomial of the basis Bj^ one has Q n (Qi(t)) = 0,Vi > n. 
Consequently, the undetermined coefficients method applied to the basis 
Bn leads always to an overdetermined system of n + 1 equations and n 
unknowns, of triangular type, which we call the fundamental system of the 
rule. We denote this system by T = .F/V",B := {-^^ = c} emphasizing that 
the system is univocally determined by the set of nodes M and its associated 
basis B. In matricial form, the fundamental system F oj = c is given by 



1 1 

01,2 















1 

h,3 
h,3 






1 



■>n-l,n 





and 



(n+l) xn 



/to 

/'I 

L 



(n+l)xl 

where //q = /(q^+i) and w = [wi, . . . ,w n ] T . The entries in F are 0jj = 



Jh 



1, . . . , (n - 1), j = i + 1, . . . , n. Since (f> (ti, 



1 and 



= (tj+i — U) 7^ 0, 1 < i < (n — 1), the matrix F has rank 

n, and so the recursive scheme ([5]) holds. 

2.1 The least-squares solution of the fundamental system 

The fundamental system ([7]) can be written in the following form 

A 



Fuj 



0...0 



(8) 



where A is the (upper triangular) submatrix of F obtained by deleting its 
last row. 

It is easy to show that the solution uj = A c is a least-squares solution V 
of the fundamental system ©. Indeed, the least-squares solution of (jSJ) is 
the solution of the normal system F T F y = F T c, and 










F T F = 


A T 












A 
0...0 



A T A, 



nx (n+l) 



(n+l)xn 
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and 










F T c = 


A T 












nx (n+l) 



Then, 



F T Fy = F T c 



HQ 



A T Ay = A 1 c. 



Ac. 



(n+l)xn 



0) 



As the system Au = c has a unique solution, this solution coincides with 
its least-squares solution. So, from ([9]) it follows that the weights's vector oj 
coincides with the least-squares solution of the fundamental system ([8]) . 

n+1 is defined as 



Recall that the p-norm of a vector x 6 

n+l 

\p ~ 



(p> i) 



The limiting case is the norm ||x||oo = maxi<i<r n+ i\ \xi\. 

The next result gives an essential property of the residual vector at the 
least-squares solution of the fundamental system F = Fj^^b- 

Proposition 2.1. Consider the interpolatory rule (pQ) and Fx = c its fun- 
damental system (jSJ). 

(a) The vector of the weights oj of the rule is the least-squares solution of 
the fundamental system. 

(b) For any p- norm, the residual vector at the least-squares solution, r(oj) = 
Flu — c, satisfies the equalities 



\\r{u)\\oc = \\r(oj)\\ p = \ij,q\, p = 1,2, 
Proof. It remains to show the equalities in (jlOp . Since 



A l.i — n 1 

r(co) 



(10) 



o 



the result follows trivially by definition of the p-norm. 



□ 



Remark 2.2. The p-norms are monotone in the sense that ||x||oo < ll x lli ^ 
1Mb < 1Mb ^ • • • ( see E] P- 43), thus the equalities observed in (fit)]) are 
not verified in general. As referred before, F is equivalent to the (2 n) x n 
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system ([2]) which arises when the undetermined coefficient system is applied 
by brute force to the standard polynomial basis < 1, x, . . . , x 2n ~ l >. So, the 
(n+ 1) equations in the fundamental system identify the relevant equations 
of the system ([2]), capturing the essential characteristics of the particular 
quadrature rule, namely its degree. Furthermore its residual vector r{u) 
satisfies the exceptional conditions given in Proposition 12.11 

2.2 The minimax solution of the fundamental system 

It is well known that for a general overdetermined linear system a minimax 
solution is not easy to compute. However, when the system is a full rank 
system of n + 1 equations in n unknowns its minimax solution can be easily 
obtained from the least-squares solution. We start by stating a result in [6] 
which relates the minimax and the least-squares solutions of this type of 
systems. 

In what follows, given a system of n + 1 equations in n unknowns F x = c, 
we denote by r(x) the residual vector r(x) = F x — c. 

Theorem 2.2. (f^, p. ^1) Let x be the least-squares solution of a system of 
n + 1 linear equations in n unknowns: rj(x) = 0, i = 1, . . . , (n + 1). Assume 
that the system is of rank n. Then, the minimax solution is the exact solution 
of the system ri(x) = Oi e, i = 1, . . . , (n + 1), where <Tj = sign(rj(x)) and 

c _ \\r{x)\\l _ i:tl\nii)\ 2 (n) 
Hr(i)||i Efflfi&l" 

We now apply Theorem l2.2l to the fundamental system F x = c given in ([8]), 
in order to obtain its minimax solution. The components of the residual 
vector at the least-squares solution u are 

r n+ i(u) = [iQ, n(u)=Q, l<i<n. 

Assuming the convention sign(0) = 1, we define the following vector of signs, 

a = [ax,.. . ,o n ,o n+1 } T = [1,. . . ,l,sign(-^ Q )] T = [1, . . . , 1, - sign(^ Q )] T . 

By Theorem 12.21 the minimax solution of the fundamental system is the 
(unique) solution of the system Ti{z) = o~i e. That is, the solution of the 
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system 



z\ +z 2 H +z n - e = fi 

01,2 Z2 H +<t>l,n z n - e = |Ul 

■-. •-. (12) 

(ftn— l,n Zn (■ — /•*?»• 1 

sign(^g) e = MQ- 



Solving for e the last equation of (fl2j) . we conclude that the required positive 
value e is simply the magnitude of the moment fj,Q, thus 

e = sign(/i Q )^ Q = |/i Q | > 0, 

as it should be by (jlip in Theorem 12,21 Substituting the value of e in the 
first n equations of (|12p we obtain the following (upper triangular) system: 

Az - \hq\v = c, where v = [1, . . . , 1] T and c = [/x , . . . ,/i„_i] T . 

The minimax solution z of the system Fx = c is then the unique solution 
of the system A z — |/zq| v = c. So, z satisfies 

A z -\hq\ v = c. (13) 

On the other hand, the least-squares solution ui of the fundamental system 
satisfies Alo = c. Subtracting this equation from f)13|) . we obtain 

A(z -u) = \hq\v<=^ At = \hq\v. 

We remark that at the minimax solution z all the residuals have the same 
magnitude e = |mq|, since by (fT3|) we have r(z) = ev. So, 

lk(l)||oo = \vq\ \\v\\oo = \hq\- 

Furthermore, by Proposition l2.U (b) it follows that HK^IIp = ImqI> f° r an y 
p-norm. Thus, ||r(io)||p = ^(z)!)^. 

The next proposition summarizes what we have just proved. 

Proposition 2.2. The fundamental system ([8]) has the following properties: 

* * 

(a) The minimax solution z is equal to z= oj + t, where u is the least-squares 
solution of the system. The vector r is the solution of the following upper 
triangular system 

At = \liq\v, withv = [1,...,1] T el", (14) 
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and hq is the principal moment of the rule. 

(b) The residual vectors at the minimax and least-squares solutions satisfy 
ImqI = ll r Hllp = 1 1^(5) | |oo, P € {oo,l,2,...}. (15) 

Example 2.1. (Simpson's rule) 

Consider the integral 1(f) = f(x)dx. The Simpson's rule is the most 
celebrated quadrature rule which uses n = 3 equally spaced nodes, say 
M = {ti,t 2 ,t 3 } = {-1,0,1}. So, 

Q 3 (/) = "1 /(-l) + L0 2 /(0) + 6^3 /(!)• 

The polynomials defining the canonical basis B are: 4>q(x) = 1, (f>x(x) = 
x + 1, 02^) = (x + 1) x, g3(x) = 02^) ( x — 1) = + 1) 35 ( x ~ 1) an d 
94(0;) = (73(2;) (x + 1) = (x + l) 2 x (x — 1). The respective moments are 
Mo = !_ 1 4>Q{x)dx = 2, m = j\4>i(x)dx = 2, fi 2 = f_ 1 ^ 2 (x)dx = 2/3, 
M3 = /_! qz{x)dx = 0, fiQ = /_! ^(x)^ = -4/15. 

Since /iQ is the first non null moment with index greater than 3, it follows 
from Theorem I2.1l -(b) that this rule has degree d = deg(Q^(f)) = 3. The 
fundamental system F uj = c for Simpson's rule is 



Fuj = 





' 1 1 1 " 








2 




1 2 




U\ 




2 




2 
_ _ 




U! 2 
. W 3 . 




2/3 
-4/15 



The weights are easily computed (recursively) by solving the triangular sub- 
system Auj = c: 

0J3 = 1/3 

oj 2 = 2-2w 3 = 4/3 

oji = 2 - (u>2 + CJ3) = 1/3. 

Thus, 

Qs(f) = ~[/(-l) + 4/(0) + /(l)]. 

Of course, by changing the domain of integration from [—1,1] to [a, b] one 
gets the traditional form of this rule: Q?,(f) = h/3 [/(a) + 4 /[(a + b)/2) + 
fib)}, where h = (b-a)/2. 
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One can easily confirm the result of Proposition ^. 11 - (a) by solving the normal 
system F T F y = F T c, that is, 



"111" 




" 2/1 " 




2 


1 2 3 




2/2 




4 


1 3 9 




. 2/3 . 




22/3 



Let us now compute the minimax solution from the least-squares solution 
co = [1/3, 4/3, 1/3] T . By Proposition 12 .21 (a), the minimax solution z of the 

fundamental system satisfies z= w+r, where r is the solution of A r = \fJLi\ v. 
So, the vector r is the solution of the system 



At 



4 
15 





" 1 1 1 " 




n 




" 4/15 " 




1 2 








4/15 




2 




. T 3 . 




4/15 



giving r = [2/15, 0, 2/15] T . Thus, the minimax solution is z= uj + r = 
[3/15, 4/3, 3/15] T . One can easily confirm the result of Proposition 12.21 (b) 
on the residual norm: 



\r(z) 



l/M| \\V\\oc 



4/15 = \\r(u) 



The Simpson's rule angle, that is the angle between its least-squares and 
minimax solutions, is 



ang(z, uj) 



43 



arccosi 



V 3V209 



180 „ r 
x ~ 7.5 

7T 



degrees. 



In paragraph 13.1.11 we pursue with other closed Newton-Cotes rules with a 
greater number of nodes. 



3 Some relevant parameters 

The properties of the least-squares and the minimax solutions of the funda- 
mental system Fj\f of an interpolatory rule Q n (f), discussed in the previous 
section, suggest the study of several features of a quadrature rule by con- 
sidering the behavior of certain parameters associated to the fundamental 
system. The first question is know what are the suitable parameters for 
previewing the accuracy of a given rule, its convergence, or even to compare 
distinct rules for the same integrand and weight functions. 

We define next a few relevant parameters associated to a fundamental sys- 
tem J-jy and we justify their relevance. In paragraph 13. II we compute these 
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parameters for a certain number of nodes using some classical quadrature 
rules. To which extend these parameters classify a quadrature rule (or dis- 
tinct rules), and the quest for other useful parameters within this least- 
squares/minimax setup, will deserve further studies. 

To emphasize the dependence on n whenever needed we index by n the 
quantities. As shown in Propositions 12. l l and 12.21 - (b) . for each n, we have 

\f i Q ) n \ = \\ r (. UJ n)\\oc = \\r(uJn)\\p = \\r(z n )\\ 00 , p>l. (16) 

These equalities reveal that the "quality" of the least-squares and minimax 
solutions is measured by the principal moment of the rule. In fact, the 
smaller is the principal moment the lesser is the magnitude of the residual 
vectors r{uj n ) and r{z n ) at the least-squares and minimax solutions respec- 
tively. 

It is well known that an interpolatory rule is convergent if and only if cj n 
satisfies the following boundness condition ([5], p. 267 or [17] . Ch. 12) 

\\u n \\i<K, n = l,2,... (17) 

This result suggests that in addition to the magnitude \fJ-Q n \ one should also 
consider the following norm parameters 

N Un = \\u> n \\i and N* = || z n ||i. (18) 

Zn 

Of course, in (]18p one can consider any other p-norm due to the equivalence 
of norms in W 1 . 

In Figure [1] we illustrate the practical relevance of the parameters N Un and 
N » for the Newton-Cotes and the Gauss-Legendre rules, with the number 
of nodes varying from n = 2 to n = 15. In the graphics of Figure Q] the 
thicker polygonal line connects the values of N UJn (represented by the dots) 
whereas the lighter line connects the values of N * . It is clear from the 

Zn 

respective graphics that for the Newton-Cotes rules the values of N LOn and 
N * do not decrease for n > 9, and so this rule cannot be convergent. We 

Zn 

shall also remark that the minimax parameter N , attains a minimum at 

Zn 

n = 9. This phenomenon is certainly related to the occurrence of negative 
weights for n > 9 in the Newton-Cotes rules ([I], p. 269), however a detailed 
study of this phenomenon is out of the scope of this work. 
If a rule Q n (f) is convergent, that is Q n (f) — > /(*) w (t) dt, then by the 
error formula © in Theorem 12. 1\ we have |/iQ„| — > 0. So, by f)16|) we get 

H^nlloo = |k(w n )||p — >0 and ||r(4)||oo — ► 0. 
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Figure 1: Comparison of the norm parameters N Wn and N * (polygonal line) 
for Newton-Cotes (left) and Gauss-Legendre rules (right), with 2 < n < 15. 



Then, in the case of a convergent rule, one should also have 

T~n = z n > 44> Z n > LO n . 

n n 

One way of measuring the closeness of oj n and z n is through what we call the 
angle of a rule Q n {f)i which we define as the angle between the minimax 
solution z n and least-squares solution co n : 

I * \ I I <z ni^n > | J / 10 n 

ang(z n ,uj n ) = arccos I — ; I , (19) 

\|| z n lb H^nlb / 

where <, > denotes the standard inner product in R n . From what we said 
before, a rule whose angle does not decreases from a certain number of nodes 
onwards cannot be a convergent rule. This is the case of Newton-Cotes rules 
(see paragraph 13. Lip . 

In Table Q] we compare the values of some parameters for Newton-Cotes 
(NC), Clenshaw-Curtis (CC), Fejer (F) and Gauss-Legendre (GL) rules, for 
n = 17 nodes and the weight function w{x) = 1. These rules are discussed 
in more detail in Section 13.11 The symbol dn denotes the degree of the 
respective rule, and the error coefficient an, taking into account ([6]), is 
defined in (120p . The computations have been carried out in a laptop using 
the Mathematica system and double precision arithmetic. 
In Table [H there is a remarkable evidence for the angle of a rule to be a 
good measure for its convergence. In fact, Table Q] shows that the angle of 
the last three rules is close to zero while the angle for the Newton-Cotes 
rule is far from zero. This means that the minimax and the least-squares 



13 



Parameters 


NC 


F 


CC 


GL 


dn 


17 


17 


17 


33 




-0.000129 


-1.07 x 10- y 


1.26 x 10~ 8 


1.80 x 10~ iU 


° 17 = A L 1 


-1.76 x 10" 20 


-1.67 x 10" 23 


1.97 x 10" 24 


6.11 x 10" 49 


ang(z, cj) 


4.55 


0.00711 


0.0380 


0.000154 



Table 1: Global comparison for NC, F, CC, GL, with n = 17 nodes. 

solutions are asymptotically aligned for the F, CC, and GL rules, whereas 
this does not happen for NC. This probably explains why in spite of its low 
degree of exactness, the Fejer and Clenshaw-Curtis rules have asymptoti- 
cally comparable accuracies to the Gaussian rules. This surprising fact has 
attracted the attention of many authors, in particular after the paper |20j of 
Lloyd N. Trefethen, where it is shown that the CC and GL rules have close 
errors for sufficiently large n. The angular values in Table [T] indicate that 
the angle of a rule deserves to be considered as (new) parameter useful to 
explain the previously referred phenomenon which holds for several rules: in 
spite of having very different degrees of exactness asymptotically the rules 
exhibit almost the same accuracy. 

For sufficiently smooth integrand functions /, Theorem 12.11 (c) gives the 
formula for the error in terms of the principal moment and of the degree of 
the rule. This justifies to consider the following parameter which we call here 
the coefficient of error. We define the coefficient of error of a rule Q n (f), of 
degree d, as 

= (W + l)!' (20) 
where \iQ n denotes the principal moment of the rule. We say that the rule 
Q n is better than the rule R n , if ctQ n < a,R n . 

Remark 3.1. Another reason for not using unbounded rules (i.e. those for 
which (fT7]) does not hold), such as the Newton-Cotes rules, is of computa- 
tional nature. When the values /(£«), 1 < i < n, have an error, say ej, it is 
easy to prove that the error EQ n in the rule satisfies 

\\e q J\ < IMIIMl, 

and so, one prefers a rule with a small norm \ \co\\. 
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For the sake of completeness, in the following proposition we relate the 
minimax and the least-squares solutions with the condition number of the 
(triangular) submatrix A of the fundamental system. 

Proposition 3.1. Letuj and z be respectively the least-squares and the min- 
imax solutions of the fundamental system (|8|) of the rule Q n (f), and I^Q its 
principal moment. Then, 



WAW II II 

\H Q \ < O, with CI = 11 111 Z (21) 

\/n 



(a) 



(b) The parameter T = - °°1 — — satisfies 

1 < r < condooiA), (22) 

where cond^A) = ||^4||oo 1 1 1 1 1 oo is the condition number of the 
matrix A for the subordinate oo-norm. 

Proof. From (j!4H in Proposition 12.21 we have 

A(z-oj) = \» q \v, (23) 
where v is the vector v = [1, . . . , 1] T £l". Then, 

* 

s/n\HQ\ < ||-4||i z Wx. 
For (b), again by ([23]) . we have r =z —oj = \{1q\ A~ 1 v, and so 

IMloo < ||^4 _1 ||oo \hq\ and |//q| < ||^4||oo [M|oo- 
The result follows from these inequalities. □ 

Remark 3.2. In the inequalities (I2ip and (I22p the influence of the nodes is 
present through the matrix A, while the influence of the functions f(x) and 
w(x) is mainly captured by the principal moment of the rule. The parameter 
r in (|22p is useful to control numerical instabilities in the computation of the 
least-squares and minimax solutions of the fundamental system, a system 
which is in general ill-conditioned due to the fact of the matrix A being in 
general almost singular for large n. These issues are however out of scope 
of the present work. 
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3.1 Worked examples 

In this section we compare some of the parameters defined in the previous 
section for the Newton-Cotes (NC), Fejer (F), Clenshaw-Curtis (CC) and 
Gauss-Legendre (GL) rules. All the computations were carried out in a 
laptop using the system Mathematica and machine precision. 

3.1.1 Newton-Cotes rules 

Previously (see Figure [Q and Table [T]) the parameters N u>n and N * were 
tested for the Newton-Cotes rules. Here we compare the angle of these rules 
as well as the norm Ht^Hoo (recall that r = u) n — z* n ), for a number of nodes 
2 < n < 17. 



70 
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40 






6 






30 


.A 




4 






20 
10 
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2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 




2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 



Figure 2: Newton-Cotes rules: on the left the angle ang(z n ,Lo), expressed in 
degrees, and on the right the distance | |t| |oo = II z n — Wn\\oo- 

The graphics in Figure [2] show that the angle between the minimax and 
least-squares solutions, ang(z n ,uj n )), has small variation for n < 8 but sud- 
denly oscillates widely. This remarkable behavior is closely related to the 
undesirable numerical properties of the Newton-Cotes rules due to the oc- 
currence of negative weights for n > 9. Moreover, from n = 9 onwards 
the distance t||oo, between the minimax solution z n and the least-squares 
solution oj, increases and so the rule cannot be convergent. 

3.1.2 Clenshaw-Curtis rules 

The Clenshaw-Curtis (CC) [7] rules are defined in [—1,1], being the points 
±1 fixed nodes. For n > 1, we consider the "practical" abscissas (see [8], p. 
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86 ) defined by 



' kir\ 

t k = cos ( — I 



k = 0, . . . , n. 



The total number of nodes is N = n + 1. Let us briefly recall our least- 
squares/ minimax approach giving the calculations for this rule when n = 3. 
First, we sort the nodes in increasing order to obtain the canonical set jV, 
that is, we take ti = — 1, t 2 = —1/2, £3 = 1/2 and £4 = 1. The respective 
triangular system is 









2 




UJ 2 




2 




UJ 3 




5/3 




UJ4 




.1/6. 



and uj 



1111 
1/2 1/2 2 
3/2 3 
3/2 



[1/9,8/9,8/9, 1/9] T its solution. So, the rule is 

CC 4 (/) = £[/(-!) + 8/ (-1/2) + 8/(1/2) + /(!)]. 



(24) 



As the principal moment is ^cc A = 1/15 ^ this rule has degree d = 3, 
like Simpson's rule. However, this rule may be considered "better" than 
Simpson's rule since its moment is smaller (fJ-cd — 1/15 < = 4/15). 




Figure 3: Clenshaw-Curtis rules: the dashed line joins the values || z Tl 
and the other line the values ||o; n ||oo- 



In Figure [3] we display the values of the oo-norm of the least-squares solution 
versus the minimax for the referred CC rules. Here we observe an interes- 
ting behavior when compared with the Newton-Cotes case: the magnitude 
of the solutions oo n and z n become closer as n increases. This suggests that 
the Clenshaw-Curtis rules are convergent, as in fact they are. 
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Figure 4: Clenshaw-Curtis rules: on the left the angle ang(z n ,Lo n ) and on 
the right H"^!^ = || z n -w n ||oo- 



The Figure S] shows values of the angle of the CC rules as well the distance 
(in the maximum norm) between the least-squares and minimax solutions. 
It is observable an enormous reduction and smoothing in the referred angle 
in comparison with the Newton-Cotes case. Here, as the value of n increases 
the angle reduces almost to zero, suggesting that both the solutions u n and 
Zfi are asymptotically coincident. 



3.1.3 Fejer rules 

The Fejer first rule [12] with n nodes uses as abscissas = cos ^^" 1 tt), k = 
1, ... ,n, which are the zeros of the Chebyshev polynomial T n (x) = cos(n 
arccos x), in [—1, 1] ([3], p. 539). For n = 3, the nodes are t\ = \/3/2, t 2 = 
and t% = V3/2. The moments are no = 1, \i\ = V^,^2 = 2/3, and 

M3 = J (t + VS/2) t (t - V3/2) dt = 

/i 4 = y {t + V3/2) 2 t(t-V3/2)dt = -1/10. 

The solution of the rule's fundamental system is [4/9, 10/9, 4/9] T . Its degree 
is obviously d = 3, like in the Simpson's rule. As the principal moment of 
this rule, |^^| = = 1/10, is smaller in modulus than the principal 
moment of the Simpson's rule, we can say that Fejer's rule is better than 
Simpson's rule. 

In Figure[5]we compare the norms ||w n ||oo arid || z ||oo, for Fejer's rules. The 
computed values suggest that Fejer's rule is convergent and the respective 
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weights are all positive. Furthermore, contrarily to NC case, as ||cj n ||oo 
— ^ II ^ 1 1 oo one infers the convergence of Fejer's rule. Both ||w n ||oo 

and 

II z n \\oo are much smaller than their Newton-Cotes counterparts, showing 
the superiority of a Fejer rule when compared with a Newton-Cotes rule 
with the same number of nodes (see also Figure [1]). 



11 13 15 17 




11 13 15 17 



Figure 6: Fejer rules: on the left the angle ang(z n ,ui n ) and on the right 



-UJr 



The Figure[6]displays the angle ang(z n , uj n ), for the Fejer rules and the values 
for the respective least-squares/minimax distance ||t||oo. Both parameters 
approach zero, as the number of nodes increases, more rapidly than the CC 
rules. 
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3.1.4 Gauss-Legendre rules 

The Gauss-Legendre rule (GL) [TU] uses as nodes the zeros of the n Legendre 
polynomial P n (x) in [— 1 , 1] (see for instance [1] , |16j ) . Gaussian rules are 
widely used in numerical quadrature in part due to their maximal degree 
d = 2n — 1 of exactness (see pT] for a survey on Gaussian rules). 

For 2 < n < 17, we compare in Figure the norms of the respective least- 
squares and minimax solutions, and in Figure [8] are displayed the respective 
angle values and distances 1 1 1 1 oo • 

It can be observed that both solutions u n 
and z n become very close for n > 5, and their distances approach zero more 
rapidly than in all the previous examples. This is of course a confirmation 
of the Gaussian rules excellence. 




o.o [ 

3 5 7 9 11 13 15 17 



Fi gure 7: Gauss-Legendre rules: the solid line represents ||wn||oo the other 




Figure 8: Gauss-Legendre rules: on the left the values of ang{z n ,uj n ) and 
on the right 11x1100 = 11 z n -w n ||oo (right). 
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